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We present a general formalism for the evaluation of time-reversal odd parton distributions. This 
formalism is applied to the evaluation of the two T-odd TMDs allowed at the twist-2 level, i.e., 
the Sivers and the Boer-Mulders functions. We have performed the calculation for two different 
models of proton structure: a non relativistic constituent quark model and the MIT bag model. 
The results obtained in both models are compared. So are the two T-odd functions. We comment 
on the fulfilment of the Burkardt sum rule as well as on the current status of T-odd TMDs in 
phenomenology. 
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Our knowledge on the hadron structure is incomplete. We, for example, know that the parton 
distribution functions describe the 1-D structure of hadrons. At leading order, the pdfs are three: 
number density, helicity and transversity. However the experimental knowledge on the latter is 
rather poor as it is a chiral-odd quantity not accessible through fully inclusive processes. Therefore, 
Semi-Inclusive DIS experiments, where a final hadron is detected, have been proposed to extend 
our knowledge on this function. This extension from fully inclusive to semi-inclusive processes 
implies a generalization of the distribution functions. That is, if one wants to study the transverse 
momentum distribution of the produced hadron, one has to account for transverse motion of quarks. 
In fact, we now know that non-perturbative effects of the intrinsic transverse momentum kj of the 
quarks inside the nucleon may induce significant hadron azimuthal asymmetries [|l|, ||]. 

It is in this context that the Sivers and the Boer-Mulders functions were defined. Transverse 
Momentum Dependent pdfs (TMDs) are the set of functions that depend on both the Bjorken vari- 
able and the intrinsic transverse momentum of the quark. 1 Their number is fixed by the number 
of scalar structures allowed by hermiticity, parity and time-reversal invariance. However, the exis- 
tence of final state interactions allows for time-reversal odd functions [pO. In effect, by relaxing this 
constraint, one defines two additional functions, namely, the Sivers and the Boer-Mulders function. 
These functions are related, respectively, to single spin and azimuthal asymmetries, and are there- 
fore important in our quest for the understanding of the proton spin. 

The Sivers function, / 17 ?(x,£r) [Qj, and the Boer-Mulders function, h i [ q {x,kj) [Qj, are for- 
mally defined, according to the Trento convention [^j, for the quark of flavor q, through the follow- 
ing expression 2 : 

f ±cj (x k T ) - -— f d ^ d2 ^ T e -K*P + S--h'h) 

f lT {x,k T )- (2?r)3 e 

*l I 5,(P5,|v7 9 (r,lr)4(oo,r)7 + ^oKO)^(0,0)|P5 v )+h.c., (1) 

z s,.=-i,i iT 

taking the proton polarized along the y axis; and 

1 V ' ; 2k x J (2tt) 3 

£ (P5 z |v7 9 (r,lr)^KD7 + 7 2 75^oKO)^(0,0)|P5 z )+h.c, (2) 

z s_=-i i i Sr 

where 5 is the spin of the target hadron. The normalization of the covariant spin vector is S 2 = — 1 , 
M is the target mass, is the quark field and S^t is the gauge link. 3 The gauge link contains 

a scaling contribution which makes the T-odd TMDs non vanishing in the Bjorken limit. 

The difference between the two functions is clearly physically transparent from Eqs. ([[]) and 
(||). The BM function counts the transversely polarized quarks, hence the Dirac operator y + y 2 /5 in 

'They also depend on the scale Q 2 , like the pdfs. 

2 a ± = (ao±a 3 )/V2. 

3 The gauge link is defined as Jzf| («j,<!;~) = ^exp (— ig J|l A + [t]^ ,l,j)dr\^ j , where g is the strong coupling 
constant. This definition holds in covariant (non singular) gauges, and in SIDIS processes, as the definition of T-odd 
TMDs is process dependent. 



2 



T-odd TMDs in Quark Models 



A. Court oy 




Figure 1: Contributions to the T-odd TMDs. 



Eq. @) in an unpolarized proton. On the other hand, the Sivers function counts the unpolarized 
quarks, hence the Dirac operator y + in Eq. ([I]), in a transversly polarized proton, i.e. the explicit 
transverse component S y in Eq. ([j]). If there were no scaling contribution of the gauge link, the two 
T-odd functions f^(x,kj) and h x q (x,kr) would be identically zero. 

1. The Interaction: MIT bag vs. NRCQM 

The general formalism for time-reversal odd functions in quark models has been presented 
in Ref. [Jsj] . This formalism implies the development into free-quark states and requires that we 
work within the Impulse Approximation. None of the quark model we have used contain gluonic 
degrees of freedom, therefore the Final State Interactions are consistently introduced through a One 
Gluon Exchange, as depicted in Fig. |l[ The expression for the Sivers and the BM functions are 
similar in any quark model: the Dirac structure dictates the spin structure of the quark operator 
defining the interaction V mm2mjm ; then the resulting spin-flavor-color matrix elements dictate the 
spin combinations allowed by such a structure and reflects also our modelling of the proton. For 
instance, in the MIT bag model [0], one gets, 4 

k x j (zn) j [zn) mu , n2 ^ m ^ m4 

(1.1) 

where no further recoil of the target is considered. The imaginary part is taken only for the Sivers 
function ; it is a consequence of re-expressing the transverse proton spin component in an helicity 
basis. The spin combinations at the quark level appear clearly from the expressions 



E T ij T kl ( PS z — ^\ b qm^ i qm 2 ^ m J :)l ^n M \ PS z 



2 



(1.2) 



4 The expressions as well as the subscripts in a [] refer to the BM function. 
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whose calculation is performed, here, assuming 51/ (6) symmetry, for q = u,d. The interaction 
V{k,qj) is here evaluated using the properly normalized fields for the quark in the bag [0] given in 
terms of the quark wave function in momentum space, which reads 

<p m {k) = iV^NRl(jf^ ) , (1-3) 
\a-kti(\k\)Xm) 

with the normalization factor N. The interaction is then 

V(k,k 3 ,q T ) mim2}tn3m = — (pi (k - q T ) fy + T m (p m2 (k)(pl (k 3 )fy + (p m (k 3 - q T ) , (1.4) 
q 

where T f W = 1 or y 2 /5 for, respectively, the Sivers and the Boer-Mulders function. The expressions 
in the NRCQM are similar to the one in the bag model. They read, with the intrinsic proton wave 
function, 

2 ME P f d 2 q T ^ 



f 1T [hl] Xq (x,k T ) = -2ig 2 ^ J 



(2k) 5 



J dk,. dk 3 (2K) 3 8(k+ -xP + )8(k l± +q ± -k x ) 
9 L L L 5 « N ^Is (*l ' i m i ' ^} ; ^3 , {m 3 ,k,P};-% 3 -ki, m n ) T-jTfi 

1 S p ijkl 

V(ki,h,q)mim 2 ,m 3 m^rS (%1 + <Z> { m 2, j, q}\ % ~ Q, {m,l, P}', -%-k\,m n ^ 

(1.5) 

The intrinsic proton wave function explicitly depends on the momenta so that there is no possible 
factorization of the spin-flavor-color matrix elements, like in Eq. ( |1.1[ ). The interaction here reads 

V(ki,k 3 ,q) mim2mm = \u nii (ki)Y + T f ^u mi (k\ +q)u m (k 3 )y + u m (k 3 -q) , (1.6) 
with u(k) the four-spinor of the free quark states. 
2. Results with an 51/(6) proton wave function 



In Refs. N8L |9[ |10fl , we have evaluated both T-odd TMDs in both the MIT bag model and a non 
relativistic Constituent Quark Model (NRCQM) using SU (6) symmetry for the proton. 

In the former case, the qualitative results are directly reflected from the calculation of the 



coefficients (1.2). In effect, through those coefficients it is possible to reconstruct what happens 



at the level of the quark helicity in a perfectly transparent way. In the latter case, one first has to 



re-express Eq. (|L5J) in terms of the formal expression of the proton state. In spectroscopic notation 



and with the Jacobi coordinates, one has 

e-( k l +k Dl 



\ 2 Sv,2)s = 3/2 3 \x)s , (2-1) 
where \%)s is the standard SU (6) vector describing the spin-flavor structure of the proton. 



4 



T-odd TMDs in Quark Models 



A. Courtoy 



Therefore, in the MIT bag model, the role of each contribution can be followed and evaluated. 
Similar conclusions can be driven from the NRCQM calculation but the expression are slightly 
more intricate due to the momentum dependence of the proton wave function. It is however worth 
noticing that the main difference between the two model's philosophy is that the interference term 
giving rise to a non-zero Sivers function arises either due the MIT bag wave function or to the 
4-spinors of the free quark states, in the CQM approach. This appears clearly in comparing the 



expressions Eqs. ( L4| , |1.6| ) for the interaction. There are only 2 spin combinations contributing to 
the Sivers function. The dominant contributions to this function comes from the spin-flipping of 
the quark interacting with the photon, i.e. the Y term in the MIT bag calculation 

fj#(x,br) = y^C 2 | ^± I lC^Y(q ± ,k T )+C+-U(q ± ,k T )] , (2.2) 

with C a normalization factor, C q a weighting spin-flavor-color factor resulting from the matrix 
elements (1.2) and where Y /U(q±,kr) include the momentum dependent part. 5 



On the other hand, there are more spin combinations for the BM function. 6 The first reason is 
that both non-flipping and double-flipping terms are important. The second reason is the sum over 
the two spin states, i.e. S z = — 1, 1. Due to the spin-flavor-color coefficients, i.e., due to the SU (6) 
symmetry assumption, the non-flipping term weights more than the double-flipping contribution. 
In effect, the latters are governed by the product of the two lower components of the bag wave 
function which encodes the most relativistic contribution arising in the MIT bag model. They turn 
out to be a few orders of magnitude smaller than the dominant ones, arising from the interference 
between the upper and lower parts of the bag wave function. This also happens if a proper non 
relativistic reduction of the gauge link, suitable for CQM calculations, is performed, justifying 
then the non relativistic approximation. 

On Fig. U we show both the Sivers and the Boer-Mulders functions for u and d quarks in both 
the CQM and MIT bag model, with the value a s .(/i^)/(47r) ~ 0.13 [jl2|]. We next explain why we 
choose such a value for the strong coupling constant. 



3. From Model Calculation to QCD 

Model calculations are important for their inputs to our current knowledge about TMDs. Most 
of the models present in the literature are relevant at scales as low as 0.1 GeV 2 . This scale is 
identified by matching the experimental value of the second moment of the unpolarized pdf with 
the result obtained within the model. It moreover represents the energy from which we assume 
perturbative QCD to be applicable. 



Nevertheless, the latter argument is controversial but, so far [13], that's the best we can do. 
The problem is a bit more tricky when one goes to transverse momentum dependent distributions, 
like we are doing here, as the TMD community hasn't come through the TMD QCD evolution's 
complexity yet. Basically, by trying to link model calculation with QCD, one faces the 2 problems: 
a "too" low initial scale and no workable evolution equations. 



5 SeeEqs. (8-9) of Ref. [§]. 
6 SeeEq. (13) of Ref. [|k]]. 
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Figure 2: Comparison of both the Sivers (red) and the Boer-Mulders (blue) functions in the NR Constituent 
Quark Model (left) and the MIT bag model (right). Small dashed curves represent the of-quark distributions 
; full curves the w-quark. 



In Refs. [|, g] we have nevertheless tried our hand by comparing our results for the first 
moment of the Sivers function with the extraction from the HERMES data by evolving, up to the 
experimental scale, our low-scale result. The trend followed by both curves is similar but it is clear 
that both results can be improved; e.g. the model calculation can be improved by evolving properly 
; the parametrization also, by, e.g., including more data. 

In a more naive spirit, we here give a comparison of the Boer-Mulders function calculated in 
both the MIT bag model and a NRCQM with a very first extraction from the unpolarized SIDIS 



data on the cos 20 asymmetry from COMPASS and HERMES Q14Q. The latter is plotted in green 
in Fig. H with no error bands as the authors of Ref. [|l4|] have estimated the theoretical errors to be 
bigger than the errors from the data. This justifies that they would consider neither the theoretical 
nor the experimental ones. A proper evolution of the model results to the experimental scale would 
push the distributions towards small x, as happens with pdfs. Also, we expect the magnitude of the 
distribution to slightly decrease when evolved to higher scales. However, for the two reasons we 
have given above, no conclusion can be driven concerning the behavior in Q 2 of the TMDs. 

4. The Burkardt Sum Rule 

There exist many model calculations as well as many variations on some models ; we refer to 
Refs. [|| ^, |K]] for bibliography. We nevertheless think it is worth comparing our results with a 
recent evaluation of the T-odd TMDs in a light-cone quark model [|TJ]. As it could be expected, the 
results obtained by overlap of Light Cone Wave Functions are consistent with ours, e.g., predictions 
similar in shape and non proportionality between the u and the d quark distributions. The latter 



point is of great interest as it is related to the Burkardt Sum Rule [15]. 

Apart from the scale input, model calculations are important in that some first principle prop- 
erties can either be checked or lead to constraints. By first principle properties we refer, in the 
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Figure 3: The Boer-Mulders function in both the NR Constituent Quark Model (red) and the MIT bag model 
(blue). Small dashed curves represent the li-quark distributions ; full curves the M-quark. The green curves 
corresponds to the parameterization at SIDIS'scales (see text). 



distribution functions framework, to sum rules. The Burkardt sum rule is related to momentum 
conservation: as there is no net transverse momentum in the proton, the sum over all partons of 



the first moment of the Sivers function should vanish Q15Q. In Refs. [g, |9|] we have shown that the 
Burkardt sum rule is fulfilled at a level of a few percent for both the MIT bag model and the CQM. 
This "discrepancy" cannot be fixed by correcting the support and it has therefore to be understood 
as a shortcoming of the two models. Moreover, this shortcoming should somehow be related to 
momentum conservation. As the bag states are not good momentum eigenstates, we do not expect 
the bag model results to fully satisfy the Burkardt Sum Rule. Also, the approximations that we 
have used in the CQM calculations result in a breaking of the momentum conservation. Namely, 
the Impulse Approximation as well as the non relativistic reduction forbids momentum to be fully 
conserved. This defect could not possibly be fixed by relaxing the non relativistic reduction or by 
applying the techniques of Ref. JT2[]. It is actually the Impulse Approximation that renders this 
recovery impossible in our approach. 

On the other hand, the LCWF fully satisfies the Burkardt sum rule. While it is fully relativis- 
tic, it also makes use of the Impulse Approximation, but, in this case, it allows for momentum 



conservation. The difference between our approach []8|] and the approach of Ref. J16Q lies in the 
dynamical form, e.g. we work in the instant form while the LCWF approach is defined in the front 
form. Focusing on this difference, we conclude that the results of using an approximation depends 
on the dynamical form. 
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